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The two-loop relations between the running gluino-quark-squark coupling, the gluino and 
the quark mass defined in dimensional regularization (DREG) and dimensional reduction 
(DRED) in the framework of SUSY-QCD are presented. Furthermore, we verify with the 
help of these relations that the three-loop /^-functions derived in the minimal subtraction 
scheme combined with DREG or DRED transform into each other. This result confirms the 
equivalence of the two schemes at the three-loop order, if applied to SUSY-QCD. 
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DRED has been introduced in Ref. [1] as a regularization scheme for supersymmetric gauge 
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■ theories which maintains supersymmetry (SUSY) and at the same time retains the elegant 

features of DREG [2], especially the gauge invariance. The essential difference between DRED 
and DREG is that the continuation from 4 to D dimensions is made by compactification. After 
dimensional reduction to D = 4 — 2e, it is only the D components of the gauge field that 
generate the actual gauge interactions. The remaining 2e components behave under gauge 
transformations as a multiplet of scalar fields, usually called e-scalars. 

As pointed out by Siegel himself [3], there are potential problems with DRED. In Ref. [4] it 
has been shown that the variation 5S of the action of a pure (no chiral matter) supersymmetric 
gauge theory is nonzero even with DRED. If 5S gives a nonzero result when inserted in a Green's 
function this creates an apparent violation of supersymmetric Ward identities. Within DREG 
this happens at one-loop order. On the other hand, within DRED all explicit calculations up to 
two- loop order have found zero for such insertions [5,6]. Recently, a mathematically consistent 
formulation of DRED [6] and rigorous methods to prove its supersymmetric properties [7] have 
been introduced. 

Another way to verify the consistency of DRED with SUSY is to study the behaviour under 
the renormalisation of the e-scalar-couplings (also called evanescent couplings) to matter and 
gauge fields. In a supersymmetric theory, they have to remain equal to the gauge coupling, 
if the renormalization scheme preserves SUSY. Explicit computations up to three-loop order 
within SUSY-QCD [8] confirmed this requirement for DRED in combination with the minimal 
subtraction scheme, i.e. the DR scheme. But, if DRED is applied to non-super symmetric 
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theories, like for example QCD, this equality is not preserved under the renormalization [9, 
10]. However, even in softly broken supersymmetric theories like the Minimal Supersymmetric 
extension of the Standard Model (MSSM), one has to worry about the e-scalars. In such theories, 
they will receive a loop- induced mass, which will also influence the renormalization of the genuine 
scalar masses. In order to decouple the e-scalar masses from the /3-functions of the genuine scalar 
masses, additional finite counterterms proportional to the e-scalar masses have to be added to 
the renormalized scalar masses. This new renormalization scheme, usually known as the DR' 
scheme, was introduced in Ref [11] to the one-loop order and extended through two- loops in 
Ref. [12]. The results presented in this letter are the same in the DR and DR' schemes, because 
we did not take into account dimensionful couplings. 

As is well known, the equality of the Yukawa couplings of gauginos to matter multiplets and 
the gauge couplings, or the equality of the quartic scalar couplings, e.g. four-squark or four- 
slepton couplings, and the gauge couplings are not preserved under renormalization if DREG is 
employed. This is a direct manifestation of the fact that DREG breaks SUSY. It means that, 
if one demands that the renormalized couplings are the same at some renormalization scale, 
then they are different at another scale. This point becomes important if we want to relate a 
given theory at one scale to the same theory at another scale. This procedure is often known 
as the running analysis and it amounts to determine the fundamental parameters of the MSSM 
solving the system of the Renormalization Group Equations (RGEs) with two types of boundary 
conditions: i) universality conditions imposed at some very high energy scale like the unification 
scale; and ii) low-energy constraints obtained from experiment. The appropriate renormalization 
scheme at each step of the running analysis is not fixed a priori. In general the SM parameters 
and cross sections are mostly given in the MS scheme [13], while the MSSM ones are usually 
given in the DR scheme. Apart from the finite shifts of the running parameters associated with 
the change of renormalization scheme, also threshold corrections, which account for the non- 
decoupling of heavy particles in mass independent schemes have to be implemented. They are 
known at the one-loop order for the complete MSSM [14], and at the two-loop orders for the 
SUSY-QCD [15,16]. 

Very recently, Refs. [17, 18] have shown that the QCD factorization theorem holds through one- 
loop order, if DRED is employed in computations of hadronic processes. They also provide 
translation rules from DRED to other regularization schemes through one-loop. However, it 
seems that the application of DRED to hadronic processes beyond one-loop becomes much more 
involved as compared to the standard procedure based on the DREG. It is thus advisable to use 
different regularization schemes for various parts of a practical computation. The consistency of 
such an approach is guaranteed by the fact that DRED and DREG are equivalent to all orders 
in perturbation theory if applied to a renormalisable theory [19]. This means that the two 
schemes are related by coupling constant redefinitions, under which the /3-functions calculated 
in one scheme transform into those computed in the other one. In the framework of QCD, 
the translation rules for the change from DREG to DRED is known up to three loops for the 
strong coupling constant and the quark masses [10,20]. In the MSSM, the one-loop relations are 
known for the gauge, Yukawa, quartic scalar couplings and for the coupling associated with the 
gaugino-chiral supermultiplet interactions, as well as for the gaugino masses [21]. The one-loop 
relation between the gauge coupling constant and the one associated with the interaction of 
the gluino and the quark-squark multiplet has also been verified by an on-shell computation in 
Ref. [22]. For the strong coupling constant even the two-loop conversion rule in SUSY-QCD is 
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known [15]. 

It is the purpose of this letter to extend the translation "dictionary" between the two schemes 
in the framework of SUSY-CQD to two-loop order. More precisely, we give in Section [2] the 
differences between the running gluino-quark-squark coupling and the running quark and gluino 
masses computed in the MS and the DR schemes at the two- loops. As a by-product result we 
reconfirm the two-loop conversion relation derived in [15]. In Section [3] we explicitly verify that 
the three-loop DR /3-functions and the fermion mass anomalous dimensions can be obtained 
from the MS results just converting all running parameters (couplings and masses) according to 
the two-loop results derived before. In Appendix [A] we discuss the one-loop renormalization of 
the four-squark couplings within the MS scheme. 

2 Two-loop conversion rules from DRED to DREG 

In this letter, we restrict the discussion to the translation rules for the running parameters of 
the SUSY-QCD. Thus we just need the SU(3) part of the MSSM Lagrangian. However, we give 
here the results valid for a general super symmetric theory based on an SU(N) gauge group, with 
one gauge supermultiplet in the adjoint representation (A), comprising the gluon and gluino, 
and Nf sets of matter multiplets in the fundamental representation (F), containing the fermions 
and their superpartners0 

2.1 Running coupling constants 

In order to compute the relations between running parameters defined in two different renor- 
malization schemes, one has to relate them to physical observables which cannot depend on the 
choice of scheme. For example, the relationship between the strong coupling constant defined 
in the MS and DR schemes can be obtained from the S-matrix amplitude of a physical process 
involving the gauge coupling computed in the two schemes. However, beyond one-loop the com- 
putation of the physical amplitudes becomes very much involved. We applied this method only 
for the computation of the two-loop effective charges of the gluon-quark-quark and gluino-quark- 
squark couplings in the DR scheme, in order to prove the equality of the corresponding couplings 
at this order in perturbation theory. We considered the simplifying case of a supersymmetric the- 
ory, i.e. massless gluino and equal-mass quarks and squarks and required the external particles 
to be on-shell. For the computation of the resulting two-loop on-shell integrals we used existing 
automated programs [23]. The effective charges computed for on-shell gluons and gluinos are not 
infrared safe, but the infrared divergences of the two charges are equal. This can be understood 
from the fact that they are proportional to the corresponding one- loop effective charges, which 
have been shown to be equal [22] , and the proportionality factors are universal quantities equal 
for gluon and gluinos in a supersymmetric theory. We found that the two effective charges are 
equal, which implies that the couplings themselves are also equal in the DR scheme through 
two-loops. The equality of the two couplings in the DR scheme has been confirmed even at the 
three-loop order in Ref. [8]. This result proves on the one hand the supersymmetric character 
of the DR scheme, and on the other hand it allows us to derive the relation between the two 

We work with Dirac fermions and complex scalar fields. 
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couplings valid in the MS scheme, as we discuss below. 

For the computation of the translation relations between the MS and DR schemes we employed 
a simpler computation method [10]. Starting from the observation that the ratio of the charge 
renormalization constants calculated using DREG or DRED is momentum and mass indepen- 
dent, one can derive them avoiding the use of the on-shell kinematics. Instead, one introduces 
physical renormalization constants, which are computed choosing a convenient kinematics for 
which the "large-momentum" or the "hard-mass" procedures can be applied, and retains the 
divergent as well as the finite pieces of the renormalization constants. Up to three loops this 
procedure is quite well established ( for a detail description of the method see Ref. [10]) and 
automated programs exist to perform such calculations [24-26]. 

Considering the physical charge of the gluon-quark-quark coupling at two-loop order we reconfirm 
the result derived in [15]. For completeness we reproduce it here 
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where a^ s = (g^ s ) 2 / (4tt) and a^ R = (g^ R ) 2 / (4tt) denote the strong coupling constant in the 
MS and DR scheme, respectively. We choose the usual normalization for the Dynkin index T F of 
the fundamental representation Tr(T a T b ) = T F 5 ab = \5 ab . Accordingly, the quadratic Casimir 
invariant for the fundamental representation is given by C F = T f Na/cI(F), where Na = N 2 — 1 
is the number of generators and d(F) = N is the dimension of the fundamental representation. 
The Casimir invariant for the adjoint representation reads Ca = N. 

Similarly, one can determine the conversion rules for the coupling constant a s = (<? s ) 2 /(47r) of 
the Yukawa interaction of the gluino and the quark-squark multiplet 



qL,ig a qL,j - qn,ig a qR,j + h.c. 



(2) 



Here g, q and q denote as usual the gluino, quark and squark fields, L and R subscripts stand for 
the left- and right-handed components of the quark and squark fields, and a and i,j are color 
indices of the adjoint and fundamental representations, respectively. 

Let us remark that we performed the calculation for a general covariant gauge and used the 
cancellation of the gauge parameter in the final results as an internal check. For the derivation 
of the two-loop formulae given above, also the one-loop relation between the gauge parameter 
defined in the DR and MS schemes is necessary. In order to properly take into account the 
Majorana character of the gluino, the rules given in [27] are applied with the help of a specially 
written PERL program [28]. 

So, for the two-loop conversion rule of the gluino-quark-squark coupling, we obtain 
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and together with Eq. (pQ) we get the relationship between and a^ s 
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As a consistency check, we will show in Section that the three-loop /3-functions of a s and d s 
computed in the MS scheme can be converted into the DR /3-function [8, 29] only by means of 
the finite shifts of the running couplings. 

2.2 Running fermion masses 

The particle masses are other fundamental parameters of the MSSM, that acquired a lot of 
attention both theoretically and phenomenologically. In this letter, we provide the two-loop 
translation relations for the fermion masses. They are functions only of the coupling constants 
and colour factors. The relations between the running masses defined in MS and DR can be 
obtained using the same requirement as for the coupling constants, that physical observables 
have to be renormalization scheme independent. 

In practice, we have employed the easier method of physical renormalization schemes as discussed 
above. So, the running quark mass defined in the MS scheme can be translated into the running 
mass in the DR scheme through 
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For the running gluino mass we get the following conversion relation 
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Again, one can verify the correctness of these relations by showing that the three-loop mass 
anomalous dimensions computed in the MS scheme can be translated into the DR ones, by em- 
ploying only the mass and coupling redefinitions given above. This point will be discussed in 
detail in the next section. 

Let us point out that the relations between the running masses defined in different renormal- 
ization schemes are free of the renormalon problems which affects the pole masses. It is thus 
advisable to use these relations in high precision calculations of the supersymmetric mass spec- 
trum. 



3 Three-loop renormalization group functions in DREG 

The renormalization group functions provide the scale variation of the parameters of a quantum 
field theory. They have been extensively studied and an impressive theoretical accuracy has 
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been achieved. In the MS scheme, the anomalous dimensions of all SM parameters are known 
up to two-loop level [30,31], while for QCD even the four-loop order results are available [32-35]. 
For a more general theory containing gauge, Yukawa and quartic scalar interactions, the gauge 
/3-function is known through three-loops [36] both in the MS and DR scheme. In the case of the 
MSSM, the three-loop anomalous dimensions for dimensionless as well as dimensionful couplings 
were derived in the DR scheme in Refs. [29,37,38]. The three-loop anomalous dimensions for 
the dimensionless couplings of SUSY-QCD were re-confirmed in Ref. [8]. 

In this section, we discuss the results for the three- loop /3-function of the gauge and gluino-quark- 
squark couplings and the three-loop mass anomalous dimensions of the quark and gluino masses 
in the framework of SUSY-QCD with MS as renormalization scheme. For such a calculation one 
can exploit that the divergent part of a logarithmically divergent integral is independent of the 
masses and external momenta. Thus the latter can be chosen in a convenient way: we set to 
zero all masses and one of the external momenta in the three-point functions paying attention to 
not introduce spurious infrared divergences. The resulting three-loop integrals can be evaluated 
with the help of existing programs [24,25]. At the three-loop order in perturbation theory, the 
use of 75 requires special care. We adopted here the prescription introduced in Ref [8]. 

Apart from the technical difficulties, related to the genuine three-loop calculation, one has to 
bare in mind that the couplings of the gluino-quark-squark and four-squark interactions are 
different from the gauge coupling even at the one-loop order, if the_MS scheme is employed. 
Since the four-squark couplings occur in the two-loop /3-function of ck^ s , one needs their one-loop 
renormalization constants for the derivation of the three-loop /3-function of d^ Is . In addition, 
for the conversion of this result into the DR scheme the one-loop translation rules from MS to 
DR of the four-squark couplings are needed. They have been known for quite some time for 
a general renormalizable theory with scalars, fermions, and gauge fields at one- and two-loop 
order [30,39]. In SUSY-QCD the tree-level four-squark interaction is given by 

4m = - £ - - «r?A) (7) 

A,B Z 

with A, B flavour indices, and a and i,j,k,l colour indices. At the tree-level, the four-squark 
couplings are equal to the gauge coupling. After renormalization in the MS scheme, one has 
to distinguish four types of quartic scalar couplings: i) the coupling of squarks with the same 
chirality and flavour gf; ,9r, h) the coupling of squarks with different chiralities but the same 
flavour g A R , g RL , iii) the coupling of squarks with the same chirality but of different flavours 
9l B i 9r B i i y ) the coupling of squarks with different chiralities and flavours g A R , g A ^ ■ Another 
subtlety which occurs beyond tree-level is that the group colour factors do not factorize, so that 
one has to keep track of various colour tensors in the computation of the one-loop renormalization 
constants. We introduce the following tensors for the quartic squark couplings 
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and the associated coupling constants 

We provide in Appendix [A] the one-loop MS /3-function for the coupling tensors retaining the 
complete colour structure dependence. The calculation in the DR scheme is significantly simpler 
since the colour tensors factorize. The resulting /3-functions of scalar couplings are equal to the 
gauge /^-function as required by SUSY. 

The translation rules for the four-squark couplings can be obtained from the finite pieces of 
the charge renormalization functions computed in the two schemes. We did the calculation for 
vanishing external momenta and regularized the infrared divergences giving a common mass to 
all particles [40]. To one-loop order they read: 

(Sf MS )o;H = (^' DR ) li;fei - ({T a , T%{T a , T b } kl + {T a , T h } u {T\ T b } kj ) ,5 = A,X = L, 

DR 

(Sj MS )y;H = (S S ^ R ) ij . k i--^ r {T a ,T*} ij {T a ,T*} k i otherwise. (10) 

Here {T a ,T b } denotes the anti-commutator of the group generators. 

The one-loop translation rules from MS to DR of the quartic scalar couplings are known for the 
case of identical flavour scalars [21]. These relations coincide with those of (Sfyij-^i couplings 
inSUSY-QCD. 



3.1 Three-loop /3-functions in DREG 

The /^-functions for the gauge and the gluino-quark-squark couplings are defined through 
Writing 

= X>^' W ' a = a s ,a s , (12) 
i=l 

where (i) stands for the loop order, we find for the gauge /3-function 
^8,(1) = 4{ _ 3CA + 2N f T F ), 
/ rMS\ 3 / Ms\ 2 ~MS 

^ S ' (2) = NH 4 (- 6C 'i + 8C '^/^ + 12C '^/ :7 »- l-^r) -^rieN f T F (C A + C F ), 
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In the expression for (3 as A ' as well as in all the other three-loop formulae quoted in this letter, 
we identify all couplings with a^ s . The inaccuracy induced in this way is of the four-loop order, 
so that the simplified formulae are enough to perform consistency checks of the two-loop trans- 
lation relations given in the previous section. In practice, we derived the formulae distinguishing 
between the various couplings, but the results are too long to be presented here. The three-loop 
results with complete dependence on different couplings can be obtained in electronic form from 
the author. 



The three-loop /3-function of the gluino-quark-squark coupling reads 
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where Nq = N F — 1 counts the number of quark/squark flavours B different from the external 
quark/squark flavour A. The additional colour factors occurring in the above results are defined 

as 

„ df c df c N 2 -A df cd df cd N(N 2 + 6) 

D3{F) = ^V^ = ^6i\T' Da{FA) = N A = 48 ' 

d**d** _ 18-6N 2 +N* 
Di{FF) ~ N A ~ 96iV^ ' (15) 

where d^ c , d a F cd , d^ cd are the fully symmetric rank three and four tensors of SU(N), as defined 
inRef. [41]. 

^MS,(2) . g g^ ven g rg |. ag a f unc tion of the quartic scalar coupling tensors. Implementing their 
explicit expressions ([8]) one gets the RHS of the second equality sign. The explicit dependence 

on the coupling tensors is needed for the computation of the three-loop result P^'^ as a 
function of the different types of couplings. As can be understood from the above formulae, 
when the quartic scalar coupling tensors occurring in the two-loop diagrams are renormalized, 
their one-loop renormalization functions are contracted with three different colour structures. 
We did the renormalization at the diagram level, employing the appropriate colour projectors. 
For the derivation of the three-loop results with all couplings set to be equal to a^ s , one can 
avoid the introduction of coupling tensors for the four-squark interaction. In this case the colour 
structures of the tree-level couplings are preserved under the renormalization to the one-loop 
order and so, their renormalization can be done as usual. However, for the conversion of the 
three-loop MS results to the DR scheme the introduction of the coupling tensors is unavoidable, 
because the colour structures do not factorize in the second equation of the translation relations 

cnj). 

Furthermore, it is a straightforward calculation to show that employing the conversion rules 
given in Eqs. CQ, ©, CED into the MS three-loop /3-functions ([13]) and (HU), one obtains the DR 
/3-function computed in Refs. [8,29]. Since the couplings a s and a s occur already at the one- 
loop order, their two-loop translation rules are necessary to convert the three-loop /3-functions 
from MS to DR. This is a strong consistency check for the translation rules we discussed in the 
previous section. 



3.2 Three-loop fermion mass anomalous dimensions in DREG 

In this section we provide the fermion (quark and gluino) mass anomalous dimensions within 
the MS scheme through three-loops. They are derived from the renormalization constants of the 
fermion masses, which can be calculated by decomposing the fermion self-energy into its vector 
and scalar parts and then computing the counterterms for the wave functions and masses. We 
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define the fermion (quark or gluino) mass anomalous dimensions as 



rriA dfi 2 

Writing their expansion in the perturbation theory like 

3 



^ dmA a na\ 
1A = -r-o-, A = q,g. (16) 



,MS.M 

8=1 

we have for the quark mass anomalous dimension 
/ aW\ ( ^MS\ 

T MS,( 2 ) = (-??CCt - \c\ + 7C'pNpTp ) 



7a > ( 17 ) 



/n,MS\ 2 fa MS \ ( rv MS \ 11 



~MS,(3) 
'9 



3 



V 47T / 

+ (UC A C F + 47C| + 48CUC F C(3) - 48C 2 F ({3))N F T F 



(18) 



where £ denotes the Riemann's zeta function with £(3) = 1.20206. For the gluino mass anoma- 
lous dimension we obtain 
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It is an easy exercise to verify that the three-loop MS mass anomalous dimensions given above 
differ from the ones computed in DR scheme [8,37] only by the finite shifts for coupling constants 
and masses discussed in Section [2J Let us point out that, for the conversion of the three-loop 
mass anomalous dimensions the two-loop relations for masses and couplings are needed. So, this 
provide us with another important consistency check of Eqs. (pQ EJ [6]) . 
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4 Conclusions 



In this letter we present the two-loop translation rules between DR and MS scheme for the 
running gluino-quark-squark coupling and for the gluino and quark masses. We also confirm 
the two-loop relation for the gauge coupling given in Ref. [15]. Furthermore, we prove that 
the three-loop /3-function of the gauge and gluino-quark-squark couplings and the anomalous 
dimensions of the quark and gluino masses calculated in the MS scheme can be converted into 
the known DR results, by means of these two-loop parameter redefinitions. This is a powerful 
consistency check of our two-loop results. 

As a by-product of our calculation, we derive the one-loop RGEs for the four-squark coupling 
in the MS scheme and their conversion rules to the DR scheme. 



A One-loop /^-functions of the four-squark couplings in MS 



As mentioned before, the four-squark couplings behave as tensors in colour space. Their one-loop 
/3-functions in the MS scheme read 
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with 5 = A, AB , and A = L, LR . 

The new colour tensors are defined as follows 

{Ht)ij;kl = {T a T%(T b T a ) kl + (T a T b ) kj (T b T a ) u , 
(Gth-M = {T a , T b } l:j {T a ,T b } kl + {T\ T b }n{T a , T b } kj 
(&L)ij;kl = (Sf; R )ij- mn (Si R ) nm]k i + (Si R )u ]mn (S^ R ) nm]k j 
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(-^-L )ij;kl — i^L )im;nz(<Sx )mj;kn + (&L )im;kn{^L )mj;nl > 
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C* C 



As can be easily verified, even if we identify the four types of quartic scalar interactions their 
one- loop (3- functions remain different. If in addition, one sets them equal to the gauge coupling 
and to the gluino-squark-quark coupling equal, i.e. if the DR scheme constraints are fulfilled, 
then the colour structures factorize. The resulting one-loop /^-functions for the scalar couplings 
are identical with the one-loop DR gauge /3-function, as required by SUSY. 
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